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In this article, we study an antithesis of the homogeneous Fermi system (jellium model) 
namely a single multielectron atom using the sea-boson approach. This system in addition 
to having a small number of electrons(finite system) is also extremely inhomogeneous. This 
means that the plane-wave basis is ill-suited to study this system, instead we use a basis 
consisting of the eigenfunctions of an electron interacting with the nucleus situated at the 
origin (ideal atom). We compute the energy levels of this system and polarizabilities and 
compare with the results of other approaches such as Hartree-Fock and pseudopotentials for 
the former and aysmptotic wavefunctions for the latter and with experiments. 



I. INTRODUCTION 

In our previous workS we showed how to compute the single-particle properties of a homogeneous Fermi 
system using the sea-boson method. We also generalised these ideas to include spin into the formalisma. 
However, in all these approaches the implicit assumption was that the plane- wave basis was the right one to 
study the system at hand. While this assumption is no doubt an excellent one for translationally invariant 
systems and perhaps even for systems that are nearly homogeneous such as a crystal, it is extremely ill-suited 
to study very inhomogeneous systems such as organic substances or to take an extreme example like the one 
to be studied in this article, a single atom. For a single atom, we have to choose a basis that properly reflects 
the highly localised nature of the distribution of electrons in real space. For this reason we find that it is 
best not to use the plane wave basis but its antithesis, namely wavefunctions that are eigenfunctions of the 
nonintcracting electrons interacting with just the nucleus. In the next section, we show how to write down 
formulas for the number conserving product of two Fermi/Bose fields in terms of the sea-bosons but this 
time in the localised (relatively speaking, relative to plane-waves) basis which are the eigenfunctions of the 
electrons interacting with the nucleus at the origin. Following this, we write down a model hamiltonian that 
describes a system consisting of electrons interacting with themelves and with a bosonic nucleus situated at 
the origin. This may be diagonalised exactly in what is now an analog of RPA in the localised basis and go 
on to compute the energy levels of this system. 

II. PARENT-BILINEAR SEA-BOSON CORRESPONDENCE 

Let us define H = + U(r), U{r) = —Z e 2 /r. Then its eigenfunctions may be written down as follows, 

H (pi(r) = €i(pi(r) (1) 
The hamiltonian for this problem may be cast in the following form, 

K= / dr^2^{fo-)(- — + U(r))ip(fa) (2) 
J Zm 

(7 

This may be written in the canonical basis as, 

K =^2e^{t(x)c{ia) (3) 

The seaboson is defined as follows, 

^(iV) = — i= c t(? ff )(!^)* e «(-.* 7 -') c ( t V); X* 1 (4) 
J n(ia) ^ ' 
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a~ a {ia ) = (5) 

here, 

(N)=Y / n°(ia) (6) 

ia 

Here n° (ia) is the ground state distribution of electrons of the ideal atom. For example in the case of helium 
h°(i — Is, a) = 1, for either spin, and n°(i = 2s, 2p, ...a) = 0. 

c 1 '{ia)c{i a ') = A^(i'a )a^ a (i'a) + Ay , {ia)at , («r) 



1 - n(ia)y 1 - ra(iV) ^ aj^ (?c)a llCTl (i'cr ) - y n(ta)y n(i' a') ^ a! v («io-i)a io .(?i<7i) (7) 

iitTi ii<Ti 

The density operator has a similar formula, 

n(Tcr) = c\w)c(i(j) = n°(ia)j^r +^24 itTi (ia)a iiai (ia) -^a^ii^a^iixcrx) (8) 

iici iifi 

also, N a — X)? n (* cr )- Furthermore, the sea-bosons are canonical bosons, 

N CT (^'). a JpOV)] = ° (9) 

[av (iV).at (/p')l =<5- 4 7 , i , (10) 
The above correspondence reproduces the following salient features of the free theory, 

(1) The definition in Eq.( ^) when plugged into the formula for the number operator Eq.( ^) gives an identity. 

(2) The dynamical four-point and six-point functions of the free theory are reproduced properly. That is if 
one computes the correlators of the fermi bilinears in Eq. ( 0) using the free hamiltonian in Eq. ( ||) , one gets 
what one expects. 

(3) The commutation rules are reproduced only in the "RPA" -sense. 

III. THE HAMILTONIAN OF THE MULTI-ELECTRON ATOM 

Let us write down the hamiltonain of the multi-electron atom (H = K + U), 

K= / df V^(fa)( + U(r)U(ra) (11) 

J 2m 

U = — J dr J dr' v(r — r) ^ ^tp^ (ra)ift {r a )ijj(r'a )tp(ra) (12) 

Here U(f) = —Z e 2 jr and v(f — r') = e 2 /\f — r |. 

K = E + fe-^H^K^) ( 13 ) 

,0",£Tl 

Similarly, the potential term U may be written as, 



U =2 



— f dr f dr' v(f— r) ^ ] ^ (ra)tp^ (r a )tp(r a )ij)(ru) (14) 



2 



U=\ £ V y/ ,J(ta)J(ra)c(fa)c(j<j) (15) 



%iV = / df / d7 *' ^ ~ r ')^K f )^(0^j(rVy(r') (16) 



A. Properties of the Ideal Atom in the Sea-boson Language 

In the previous section E$ was the ground state energy. It is computed as follows. Let us assume we have 
N electrons in our system. This means that we have our eigenstates filled upto i max starting from i = 
B = {i = 0, 1, 2, imax}- each i has at most two electrons with spin up or down. Let us consider an even 
number of electrons, 

N = 2 x (i max + 1) (17) 

E = 2e + 2e 1 + ... + 2e imax (18) 
The ground state of this system is annhilated by the sea-bosons, 

aiAh')\G}=0 (19) 
also the sea-boson has the propety that (a postulate if you like), 

a xXi° ) = 0; f° re i <a ( 2 °) 

This allows us to construct the first excited state of this system as 

\Ei) = aU imax „((j = i max + l)a')\G) (21) 

where the spins er, a can be anything. It says that the first excited state is degenerate. That is, you can 
take a spin up or down electron situated at i — i max and elevate it to i = i max + 1 and that electron can be 
of either spin as well. Hence the four-fold degeneracy. The energy of the first excited state is given by, 

E l = 2e + 2ei + ... + 2e Wx + (e Wx+ i - e w J (22) 

Having convinced ourselves of the correctness of the formalism as regards the ideal atom we proceed to study 
the full problem. 

IV. THE NONIDEAL ATOM 

In order to compute the energy levels of the nonideal atom, we proceed as follows. First let us rewrite the 
interaction term in the sea-boson language. 

U=\ ]T V^jpJfacfaJtfa'Wa^-NVo (23) 



^o4E%K ( 24 ) 
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(25) 



Let us now diagonalise this hamiltonian. For this we have to introduce the dressed sea-bosons. 

di*= E [ d i" a l r (j tT ')\ a l,(3< T ')- E i d ^' a iaG a ')H a (3 a ) ( 26 ) 



t,J ,(7,(7 I, J ,(J,<7 

Let the diagonaliscd hamiltonian be written as, 



h = y,o>iA„*i° ( 27 ) 

la 

ui a d Ia = E (ej-ei)[dia,ai )T 0a')]a^0a')+ E (ej - ei)[d Ia , a^tfo )]al^a ) 



+ \ E V ij/r ^(Mdi„ a Ia (ja)]+A la (ia)[d Ia , a^fr)]]^, (f*)^^, (jV)+a| v (iVjA^, (,V)] 

+ ^ E \?f[ A iM°) a iM°) +A- a {m)aX a {ia)][K, a ,{j'a )[d Ia ,a^,{j' a )} + [d 7a , a^, ( t V)]A jV (zV)] 

(28) 

w 7(7 [d 7CT , 4 ct (J<t')] = (e 3 - - e;)[d 7(7 , at a (ja')} + E [A7 CT (.yV)[d 7CT , a^JJ'a)] + A^(iV)[d 7(7 , al Ji ! a)]]A 7a (ja)S a y 



* j 



(29) 



w 7(J [d 7(J ,a 3V (icr)] = ( ej - ei)[d J(T ,a^/(io-)] - E % t y [^(jV)^, a.-^jV)] + A^(i'o-)[d Jff , a 1 -, Ji' a)]]k^{ia)5^ { 

.7 r7 

(30) 

[di„al(ja')} = A ^ a)S 2"' G a (ll) (31) 

%a Via - e j + e i 



ldia,a. a ,(i*)} = - jaK °' G a (i,j) (32) 

Via - Cj + 



The eigenvalues w 7(T that correspond to energies of the excited states measured with respect to the ground 
states are given by the solution to the following eigenvalue problem, 

» j 

here we must impose the additional constraint, 
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G a (i,j) = 0; ejK* (34) 

in order that excited states have higher energies than the ground state. Further since we have constraint 
[diu, d\ a ] = 1, the matrix G a may itself be computed by the following equation, 

Also, 

d J<T = E^' at CT (jc7)]a j(7 (jcr) - ^[d /(T , a Jff (icr)]a] (T (icr) (36) 

The full hamiltonian may then be rewritten more transparently as, 

H = E 'o + E - c 0(4 T tf i)a ?T^ + 4^ ^(J T)) 



+ E w /t4t^t + E w /i4i d n ( 37 ) 

From the above form of the full hamiltonian it is clear that having obtained the eigenvalues u)j a the next 
step would be to arrange the various excited states according to whether or not (ej — ej) > ujj a There is also 
the issue of whether we are allowed to excite any number of sea-bosons all with the same label la (since 
they are bosons). The answer is no, since they have to describe particle-hole excitations of the atom and 
the elemetary particles are electrons. In order to answer this question we have to follow closely the section 
where we showed how to study the ideal atom in the sea-boson language. There however we conveniently 
glossed over this thorny issue(sort of). But it is hoped that the reader appreciates the spirit of the argument 
and will not examine this flaw under a microscope. The whole formalism has to be taken with a grain of 
salt (a tall glass of milk, a jar of pickles and a can of diet coke !). 

The ground state energy of the system may be obtained by merely computing the expectation value of 
the kinetic energy and the potential energy separately and adding them together. 

e'o = E + E( £ j - er)«<4 T tf T)a ^ t)> + <4itf u» 



+ 2 E £%?7( A &Cfc)^0^)<a? ff (^ 

-n -at -!i (t ^ CT 



+ AjJa)Aj a (fa)(a^(ia)a ila (fa)) + AjJa)Ar, a (?a)(a^a)a^(?a)}) (N)V (38) 

Since, 

a ?<r(J Cr ) = E^OM' 4a\ d I<r ~ EN^M d I*] d L ( 39 ) 



= EMM d l]dl ~ YMi^ d ^ ( 40 ) 

From these we may deduce, 

( a ? CT (/ CT )v ff (/(T)> = - Efo^)' 4jK7 ct (/ ct ). d i°\ ( 4l ) 
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^fatyjf*)) = - $>^(M 4JM^),<^] ( 42 ) 

i cr ^ 

(4 o .(Jo-)a.7 CT (j'(r)) = ^2[a^{j(r), d Ia ][a^(f a), d Ia ] (43) 
I 

(^Aja)a^ (fa)) = £ [a |(T ( ja) , d+ J [a 7 (/a), 4] (44) 
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The rest are computational details and will be done the next time we replace this preprint on the Los Alamos 
Archive. 
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